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Tj- I Abstract 
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In the framework of the Lindblad theory for open quantum systems we deter- 
(— I I mine the degree of quantum decoherence and classical correlations of a harmonic 

Qh' oscillator interacting with a thermal bath. The transition from quantum to clas- 

sical behaviour of the considered system is analyzed and it is shown that the 
classicality takes place during a finite interval of time. We calculate also the 
^ ' decoherence time and show that it has the same scale as the time after which 

, thermal fluctuations become comparable with quantum fluctuations. 
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^ : 1 Introduction 
(N ■ 

^ I The transition from quantum to classical physics and classicality of quantum systems 
! 

O ■ continue to be among the most interesting problems in many fields of physics, for 

^ ■ both conceptual and experimental reasons [1, 2, 3]. Two conditions are essential for 

^ ■ the classicality of a quantum system [4]: a) quantum decoherence (QD), that means 
I the irreversible, uncontrollable and persistent formation of a quantum correlation (en- 
tanglement) of the system with its environment [5], expressed by the damping of the 
coherences present in the quantum state of the system, when the off-diagonal elements 
of the density matrix of the system decay below a certain level, so that this density 
matrix becomes approximately diagonal and b) classical correlations (CC), expressed 
by the fact that the Wigner function of the quantum system has a peak which follows 
the classical equations of motion in phase space with a good degree of approximation, 
that is the quantum state becomes peaked along a classical trajectory. The necessity 
and sufficiency of both QD and CC as conditions of classicality are still a subject of 
debate. Both these conditions do not have an universal character, so that they are 
not necessary for all physical models. An important role in this discussion plays the 
temperature of the environment and therefore it is worth to take into account the 



differences between low and high temperature regimes. For example, purely classical 
systems at very high temperatures are described by a classical Fokker-Planck equa- 
tion which does not follow any trajectory in phase space (for very small kinetic energy, 
compared to the thermal energy, when the probability distribution becomes essentially 
independent of momentum), so that in this case CC are not necessary. Likewise, one 
can have a classical behaviour if the coherences are negligible, without having strong 
CC (for example, in the case of a classical gas at finite temperature) and the lack 
of strong correlations between the coordinate and its canonical momentum does not 
necessarily mean that the system is quantum. On the other hand, the condition of CC 
is not sufficient for a system to become classical - although the Wigner function can 
show a sharp correlation in phase space, the quantum coherence never vanishes for a 
closed system which has a unitary evolution. Likewise, in the low temperature quan- 
tum regime one can observe strong CC. For example, in the case of a purely damped 
quantum harmonic oscillator (at zero temperature), the initial coherent states remain 
coherent and perfectly follow classical trajectories of a damped oscillator, but CC are 
not sufficient for classicahty. 

In the last two decades it has became more and more clear that the classicahty 
is an emergent property of open quantum systems, since both main features of this 
process - QD and CC - strongly depend on the interaction between the system and 
its external environment [3, 6, 7]. The main purpose of this work is to study QD and 
CC for a harmonic oscillator interacting with an environment in the framework of the 
Lindblad theory for open quantum systems. We determine the degree of QD and CC 
and the possibility of simultaneous realization of QD and CC for a system consisting 
of a harmonic oscillator in a thermal bath. It is found that the system manifests a QD 
which increases with time and temperature, whereas CC are less and less strong with 
increasing time and temperature. 

2 Lindblad master equation for the harmonic os- 
cillator in coordinate and Wigner representation 

The irreversible time evolution of an open system is described by the following general 
quantum Markovian master equation for the density operator p{t) [8]: 



dpjt) 
dt 
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H is the Hamiltonian of the system and Vj, Vj are operators on the Hilbert space of H, 
which model the environment. In order to obtain, for the damped quantum harmonic 
oscillator, equations of motion as close as possible to the classical ones, the two possible 
operators Vi and V2 are taken as linear polynomials in coordinate q and momentum p 
[9, 10] and the harmonic oscillator Hamiltonian H is chosen of the general quadratic 
form 

H = Ho + ^{qp + pq), H,^^p' + ^q'. (2) 
z Am 2 

With these choices the master equation (1) takes the following form [9, 10]: 

dp X % % 

= -^[Ho,p\ - ^(A + /i)[g,pp + pp] + ^(A-/x)[p,pg + gp] 

^''-[Q: hp]] - ^b, b,p]] + ^(b, b,p]] + b, k,p]])- (3) 



^2 ......... U',U',ru . 



In the particular case when the asymptotic state is a Gibbs state Pg(oo) = e~ M"/Tre~ fc^ . 
the quantum diffusion coefficients Dpp, Dgg, Dpg and the dissipation constant A satisfy 
the relations [9, 10] 

^ A + u. ^ huj ^ X — II h hu , . 

J^vv — Aima;coth — mr, JJoa — coth —r-=, Dj,a — 0, (4) 

2 2kT' 2 mo; 2kT' pi ' ^ ^ 

where T is the temperature of the thermal bath. 

In the Markovian regime the harmonic oscillator master equation which satisfies 
the complete positivity condition cannot satisfy simultaneously the translational invari- 
ance and the detailed balance (which assures an asymptotic approach to the canonical 
thermal equilibrium state). The necessary and sufficient condition for translational 
invariance is X = fi [9, 10]. In this case the equations of motion for the expectation 
values of coordinate and momentum are exactly the same as the classical ones. If 
A ^ /i, then we violate translational invariance, but we keep the canonical equilibrium 
state. 

The asymptotic values (Tqg(oo), (Tpp(oo), o'pg(oo) do not depend on the initial values 
<7qq{0), o'pp(O), crpg(O) and in the case of a thermal bath with coefficients (4), they reduce 
to [9, 10] 

Coo(oo) = coth——, cTjjo ( oo ) = coth——, (Jdo(oo) = 0. (5) 

In the following, we consider a general temperature T, but we should stress that 
the Lindblad theory is obtained in the Markov approximation, which holds for high 
temperatures of the environment. At the same time, the semigroup dynamics of the 
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density operator which must hold for a quantum Markovian process is vahd only for 
the weak-coupling regime, with the damping A obeying the inequality A <^ a;. 
We consider a harmonic oscillator with an initial Gaussian wave function 

where CTqq{0) is the initial spread, (Tpg(O) the initial covariance, and o"g(0) and o"p(0) are 
the initial averaged position and momentum of the wave packet. The initial state (6) 
represents a correlated coherent state (squeezed coherent state) [11] with the variances 
and covariance of coordinate and momentum 

flS TlTflOJ f} T 

Here, 5 is the squeezing parameter which measures the spread in the initial Gaussian 
packet and r, with |r| < 1 is the correlation coefficient at time t — 0. The initial 
values (7) correspond to a minimum uncertainty state, since they fulfil the generalized 
uncertainty relation 

<^..(o)Mo)-s\(o) = f- 

For 6=1 and r = the correlated coherent state becomes a Glauber coherent state. 
From Eq. (3) we derive the evolution equation in coordinate representation: 

dp ih d"^ imuj'^ 2 /2\ 



For the case of a thermal bath with coefficients (4) the Wigner distribution function 
t) satisfies the following Fokker-Planck-type equation: 

dW V dW o dW d d 

+ 9 ^ ^ + P)^^ + ^TTX ■ 10 

2 2A;r ap^ ma; ag^ 

The first two terms on the right-hand side of both these equations generate a purely 
unitary evolution. They give the usual Liouvillian evolution. The third and forth 
terms are the dissipative terms and have a damping effect (exchange of energy with 
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environment). The last two are noise (diffusive) terms and produce fluctuation effects 
in the evolution of the system. They promote diffusion in momentum p and coordinate 
q and generate decoherence in coordinate and momentum, respectively. In the high 
temperature limit, quantum Fokker-Planck equation (10) with coefficients (4) becomes 
classical Kramers equation (Dpp — > 2mXkT for X — ii). 

The density matrix solution of Eq. (9) has the general form of Gaussian density 
matrices 

-it) _ ^ l^(i±^ _ ^^(,))(, _ + ' it){q - q% (11) 



where (T{t) = agq{t)app{t) — (Jpgit) is the Schrodinger generalized uncertainty function 
[12] (aqq and app denote the dispersion (variance) of the coordinate and momentum, 
respectively, and apg denotes the correlation (covariance) of the coordinate and mo- 
mentum) . 

For an initial Gaussian Wigner function (corresponding to a correlated coherent 
state (6)) the solution of Eq. (10) is 

1 1 

W{q,p,t) = j==exp{-——[app{t){q - ag{t)y + agg{t){p - ap{t)y 



-2apg{t){q-agmp-apm}. (12) 
In the case of a thermal bath we obtain the following steady state solution for 

2kT^ 



t — > oo (we denote e = ^^)'- 



In the long time limit we have also 

W^{q,p) = -———exp{-———[rmjq^ + ^]}. (14) 
nh coth e h coth e muj 

Stationary solutions to the evolution equations obtained in the long time limit are 

possible as a result of a balance between the wave packet spreading induced by the 

Hamiltonian and the locahzing effect of the Lindblad operators. 

3 Quantum decoherence and classical correlations 

As we already stated, one considers that two conditions have to be satisfied in order 
that a system could be considered as classical. The first condition requires that the 
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system should be in one of relatively permanent states - states that are least affected 
by the interaction of the system with the environment - and the interference between 
different states should be negligible. This imphes the destruction of off-diagonal ele- 
ments representing coherences between quantum states in the density matrix, which is 
the QD phenomenon. The loss of coherence can be achieved by introducing an interac- 
tion between the system and environment: an initial pure state with a density matrix 
which contains nonzero off-diagonal terms can non-unitarily evolve into a final mixed 
state with a diagonal density matrix during the interaction with the environment, like 
in classical statistical mechanics. 

The second condition requires that the system should have, with a good approx- 
imation, an evolution according to classical laws. This implies that the Wigner dis- 
tribution function has a peak along a classical trajectory, that means there exist CC 
between the canonical variables of coordinate and momentum. Of course, the correla- 
tion between the canonical variables, necessary to obtain a classical limit, should not 
violate Heisenberg uncertainty principle, i.e. the position and momentum should take 
reasonably sharp values, to a degree in concordance with the uncertainty principle. 
This is possible, because the density matrix does not diagonalize exactly in position, 
but with a non-zero width. 

Using new variables E = + and A — q — q', the density matrix (11) can 
be rewritten as 

p(E, A, t) = ./^exp[-aE2 - + i/3EA + 2aaJt)E + - /?a„(t))A - aal{t)],{15) 

with the abbreviations 

and the Wigner transform of the density matrix (15) is 

'^<''''' = ^^e.p{- '"^^fa-"''^»^'"-"-W"' - ate - a,m- (IT) 

a) Degree of quantum decoherence (QD) 

The representation-independent measure of the degree of QD [4] is given by the 
ratio of the dispersion 1/^^27 of the off-diagonal element p(0, A,t) to the dispersion 
\j2J01 of the diagonal element p(E, 0, t) : 

^QD = 2' 
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which in our case gives 



Sgnit) = — 1=. (19) 



The finite temperature Schrodinger generahzed uncertainty function, calculated 
in Ref. [12], has the expression 

a{t) = ^{e-^^*[l -{S + YTT^-^) coth e + coth^ e] 

+e-2^* coth e ^ + — ^ - 2 coth e ^^—^ 

1 ,itsin(2m) 2rita;(l-cos(2Qt)), ,2 , 

In the limit of long times Eq. (20) yields 

(t(oo) = — coth^e, (21) 



so that we obtain 



huj 

2kf' 



5Qi5(oo) = tanh^, (22) 



which for high T becomes 



Sqd{<x^) = (23) 



We see that 6qd decreases, and therefore QD increases, with temperature, i.e. the 
density matrix becomes more and more diagonal at higher T and the contributions of 
the off-diagonal elements get smaller and smaller. At the same time the degree of purity 
decreases and the degree of mixedness increases with T. 6qd < 1 for T 7^ 0, while for 
T = the asymptotic (final) state is pure and 6qd reaches its initial maximum value 
1- ^QD = when the quantum coherence is completely lost. So, when 6qd = 1 there 
is no QD and only if 6qd < 1, there is a significant degree of QD, when the magnitude 
of the elements of the density matrix in the position basis are peaked preferentially 
along the diagonal q = q'. When Sqo ^ 1, we have a strong QD. 

b) Degree of classical correlations (CC) 

In defining the degree of CC, the form of the Wigner function is essential, but 
not its position around crq{t) and ap{t). Consequently, for simplicity we consider zero 



7 



values for the initial expectations values of the coordinate and momentum and the 
expression (17) of the Wigner function becomes 

" r {^1^1 -Pf „,„2i 



As a measure of the degree of CC we take the relative sharpness of this peak in the 
phase space determined from the dispersion Ti^JlFj in p in Eq. (24) and the magnitude 
of the average of p (po — ft-Pq) [4] : 



Sec = (26) 



where we identified q as the dispersion l/\'2a of q. 6cc is a good measure of the 
"squeezing" of the Wigner function in phase space [4]: in the state (24), more "squeezed" 
is the Wigner function, more strongly established are CC. 
For our case, we obtain 



5ccit) = ^LAl. (26) 

where a{t) is given by Eq. (20) and (Jpq{t) can be calculated using formulas given in 
Refs. [9, 10]: 

= ^e-2^*{[/^^(2cothe - S - ^^3^^) - cos(21]t) 
+u;n{6 - YTT—?^) sin(2Qi) + fjLu;{6 + ^ - 2cothe) + ^JiL=}. (27) 



5{l-r^)' ' ' ^ ' 5(l-r2) Vr 



r 



When Sec is of order of unity, we have a significant degree of classical correlations. The 
condition of strong CC is 5cc <^ 1; which assures a very sharp peak in phase space. 
Since (7^^(00) = 0, in the case of an asymptotic Gibbs state, we get 5cci^) ^ 00, so 
that our expression shows no CC at i — > 00. 

c) Discussion with Gaussian density matrix and Wigner function 
If the initial wave function is Gaussian, then the density matrix (11) and the 
Wigner function (12) remain Gaussian for all times (with time-dependent parameters 
which determine their amplitude and spread) and centered along the trajectory given 
by the solutions of the dissipative equations of motion. This trajectory is exactly 
classical ior X — /i and only approximately classical for not large X — /i. 

The degree of QD has an evolution which shows that in general QD increases 
with time and temperature. The degree of CC has a more complicated evolution. 
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but the general tendency is that CC are less and less strong with increasing time and 
temperature. 5qd < 1 and 5cc is of the order of unity for a long enough interval of 
time, so that we can say that the considered system interacting with the thermal bath 
manifests both QD and CC and a true quantum to classical transition takes place. 
Dissipation promotes quantum coherences, whereas fluctuation (diffusion) reduces co- 
herences and promotes QD. The balance of dissipation and fluctuation determines the 
final equilibrium value of 5qd- The quantum system starts as a pure state, with a 
Wigner function well locahzed in phase space (Gaussian form). This state evolves 
approximately following the classical trajectory (Liouville flow) in phase space and 
becomes a quantum mixed state during the irreversible process of QD. 
d) Decoherence time 

In order to obtain the decoherence time, we consider the coefficient 7 (16), which 
measures the contribution of non-diagonal terms in the density matrix (15). For short 
times (At ^ 1, <S 1), we have: 

2 2 

TflOJ T T UJT 

From here we obtain that quantum coherences in the density matrix decay exponen- 
tially at a rate given by 

2 2 

T T UJT 

' [^^^ + S{1^)^ '"'^ W^)^ '"'^ 6 - A - - ^-^] (29) 
and then the decoherence time scale is 

^'"^ ^ 2[A(<5 + ^) coth e + ^i{5- ^) coth e - A - - ^] ' 

The decoherence time depends on the temperature T and the coupling A (dissipation 
coefficient) between the system and environment (through the diffusion coefficient 
Dpp), on the squeezing parameter 6 that measures the spread in the initial Gaussian 
packet and on the initial correlation coefficient r. We notice that the decoherence time 
is decreasing with increasing dissipation, temperature and squeezing. 
For r = we obtain: 

^'"^ ^ 2(A + /.)(5cothe-l) ^^^^ 
and at temperature T — (then we have to take /i — 0), this becomes 

_ 1 
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We see that when the initial state is the usual coherent state (5—1), then the deco- 
herence time tends to infinity. This corresponds to the fact that for T = and S — 1 
the coefficient 7 is constant in time, so that the decoherence process does not occur in 
this case. 

At high temperature, expression (30) becomes (we denote r = ^) 
_ 1 
''^'^ " 2[A(5+^)+M^-^)]t- 
If, in addition r = 0, then we obtain 

In Ref. [12] we determined the time td when thermal fiuctuations become com- 
parable with quantum fiuctuations: 

" 2[A((5 + 5(1^) coth e + fi{5- j^^) coth e - 2A] " ^^^^ 

At high temperature, expression (35) becomes 

_ 1 

" 2r[A(5 + ^)+M5-^)]- ^^^^ 

As expected, the decoherence time tdeco has the same scale as the time td after which 
thermal fluctuations become comparable with quantum fluctuations. 

We can assert that in the considered case classicality is a temporary phenomenon, 
which takes place only at some stages of the dynamical evolution, during a deflnite 
interval of time [13]. Due to the dissipativc nature of evolution, the approximately 
deterministic evolution is no more valid for very large times, when the localization 
of the system is affected by the spreading of the wave packet and of the Wigner 
distribution function. 
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